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Abstract 

The work proposes and studies a model for one-dimensional spatially extended systems, which 
involve nonlocal interactions and finite propagation speed. It shows that the general reaction- 
diffusion equation, the Swift-Hohenberg equation and the general Kuramoto-Sivashinsky equation 
represent special cases of the proposed model for limited spatial interaction ranges and for infinite 
propagation speeds. Moreover, the Swift-Hohenberg equation is derived from a general energy 
functional. After a detailed validity study on the generalization conditions, the three equations are 
extended to involve finite propagation speeds. Moreover linear stability studies of the extended 
equations reveal critical propagation speeds and novel types of instabilities in all three equations. 
In addition, an extended diffusion equation is derived and studied in some detail with respect to 
finite propagation speeds. The extended model allows for the explanation of recent experimental 
results on non-Fourier heat conduction in non-homogeneous material. 
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I. INTRODUCTION 



The propagation of activity in spatially extended systems has attracted much attention 
in the last centuries, starting from celestial mechanics in the sixteenth century to the study 
in complex physical, chemical or biological systems in the last decades. In such systems 
the communication between the system subunits plays a decisive role. These subunits may 
be coupled locally to their next neighbors or may reveal longer-distance connections. In 
the latter case, the finite speed of interactions, i.e. the finite propagation speed, may yield 
temporal delays which affect the space-time dynamics of the system. Such effects have been 
found in several systems, e.g. in neural networks 0, 0, [5, 0], quantum devices 



porous 



101 ]. hot plasma [11] 



and non-homogeneous media 18, 9fl, in metals irradiated by laser pulses 

nn 

and electronic circuits [12|, [13| as the internet or computational networks. 
In order to describe spatial systems mathematically, partial differential equations (PDE) 
have been applied widely. However, most PDE models neglect effects caused by finite prop- 
agation speeds in the system. To consider these effects in PDE models, various approaches 
have been followed as, e.g. the introduction of additional temporal derivatives [ijj] or tempo- 
ral constant delays [l5(. However, propagation delays depend on the distance between two 
locations and thus represent space-dependent delays. In contrast to the PDE formulation, 
a natural way to consider these delays are integral-differential equations (IDE) which sum 
up all activity in a spatial domain and easily consider propagation delay, see e.g . [ltl \v\ . 
Moreover, the strong connection between PDE and IDE models is well-known [181 ] . The 
present work proposes a novel IDE model involving finite propagation delays which gen- 
eralizes well-studied PDE models in one spatial dimension, namely the reaction-diffusion 
models, the Swift-Hohenberg model and the Kuramoto-Sivashinsky model. Moreover re- 
calling the mathematical description of neuronal populations by IDEs, the proposed model 
represents a generic pattern forming model with a vast range of possible applications. 
In addition to the generalization described above, the novel IDE model allows for the exten- 
sion of PDE models in order to consider propagation delay effects. The present work shows 
that the results obtained in IDE models involving finite propagation speeds can be applied 
easily to PDE models. Hence this work extends previous studies 19|, l20f] as it clearly distin- 
guishes local dynamics from nonlocal spatial interactions and elicits the strong connection 
between well-known PDE-models of pattern forming systems and the IDE models in neural 
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field theory 
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23| . We derive novel extensions of the diffusion equation, the Swift- 



Hohenberg equation (SHE) and the Kuramoto-Sivashinsky equation now involving finite 
propagation speeds. These extensions are important if the speed of the propagating activity 
in the system approaches the propagation speed of the system. Such ultra-fast phenomena 
have been observed experimentally in solids 



24, 



surfaces 
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251 ]. in plasma [ll|] and on solid and fluid 



27]. 



The present work is structured as follows. Section [TT] introduces the integral-differential 
equation model and elicits its formal relation to partial differential equations. Then the 
subsequent section shows the reduction procedure to the PDE models and investigates the 
validity of that reduction. Eventually section IIVI studies the linear behavior of three specific 
models and examines their dependence on the finite propagation speed, while the last section 
closes the work by a brief discussion of the obtained results. 



II. THE GENERAL MODEL 



The present work considers the general evolution equation 

TV{x,t) = h[V(x,t)] 

+ / dyK(x-y)f[V(y,t-T(x-y))} 

+ / dyL{x-y)g[V(y,t-r{x-y))) (1) 
Jn 

with the scalar field variable V(x, t) and the spatial domain Q which is assumed being large 
but finite. The temporal operator T = T(d/dt) represents the temporal linear dynamics 
of an uncoupled element at spatial location x. For instance, a damped oscillator would be 
modeled as TV = d 2 V/ dt 2 +dV/ dt+V. In the following, we specify the operator by T(0) = 1. 
In addition, the self-interaction h[V(x,t)] may represent an additional nonlinear driving. 
Considering the previous example, the nonlinear damped oscillators may be modeled by 
TV = V — sin(y), i.e. h[V) = V — sin(V). Moreover, to extent the system by spatial 
coupling, K(x — y) and L(x — y) in Eq. represent the different coupling functions between 
elements at spatial locations x and y and the two coupling functions belong to two different 
spatial interaction types. For instance, in neural nets neurons may be coupled by excitation 
and by inhibition, which is modeled typically by an excitatory net and an inhibitory net 
with K(x — y) > and L(x — y) < 0, respectively. In the next section, this distinction 
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turns out out to be very important in the context of the Kuramoto-Sivashinsky equation. In 
addition, the functionals /[V] and g[V] allow for various (nonlinear) interaction types of the 
corresponding spatial interactions. Finally the propagation delay t(x — y) = \x — y\/c with 
the propagation speed c takes into account the finite time it takes the signal to propagate 
from one spatial location x to another location y. 

We point out that Eq. (JTJ is similar to previous IDE models of neuronal populations [28|, [29] , 



which can be derived from basic properties in neural tissue [19J. In addition, it extends the 
previous IDE models by considering local interactions and an additional spatial interaction 
type and thus renders the model more realistic for physical systems. As we show in the 
subsequent section, this extension proves to be very powerful. 

Finally recall the relation of IDEs to PDEs. Some previous studies remarked that IDEs 
generalize PDEs 



18 



191 ] by the identity 
f 00 dyK(x - y)S[V(y)] = f>l) ^^ (x)] (2) 

J -°° n=0 ° X 

with a nonlinear functional S and the kernel moments K n = J d7]K(i])r] n /n\ Vn G N . The 
expansion (j2J) represents an order expansion of spatial interactions whose order n indicates 
the spatial interaction range. Hence K n represent the contribution of spatial interactions 
of order n. Typically local or short-range PDE models involve spatial derivatives of up 
to second order, while PDE models with spatial derivatives of up to fourth order involve 
non-local or long-range interactions. The expansion (T5]) extends this classification to spatial 
derivatives of arbitrary order. 



III. REDUCTION TO SPECIFIC MODELS 

This section shows that Eq. ([1]) generalizes react ion- diffusion systems, the SHE and the 
Kuramoto-Sivashinsky equation in one- dimensional spatial systems. Since these equations 
assume an infinite propagation speed in the system, we set t(x — y) =0. 



A. Reaction-diffusion equations 

First let us focus on reaction-diffusion models, which allow for the mathematical descrip- 
tion of diverse spatio-temporal phenomena. We mention the heat propagation or phase 
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transitions in solids 



in biological systems 



31| , chemical reactions in fluids [32j, |33j and the pattern formation 
18| . Reaction-diffusion equations are obtained from (PQ) by choosing 
T = d/dt, g[V] = and f[V] = j\V. Further, the kernel K(x) is symmetric, i.e. K 2n +i = 0, 
and exhibits short-ranged spatial interactions with K n — > for n > 2. Then the application 
of Eq. © yields 



dV(x,t) 
di 



rj2 

h[V(x,t)]+D—V(x,t) 



(3) 



with h(V) = h[V] + K fiV and D = K 2 fi- This reaction-diffusion equation [33[] accounts for 
(non-)linear self-interactions represented by h[V] and considers local diffusive interactions 
with diffusion constant D. Hence, reaction-diffusion systems neglect spatial interactions of 
order n > 4. 



B. Swift-Hohenberg equation 



In their famous work on fluctuations near the onset of the Rayleigh-Benard convection [34, 
351 ]. Swift and Hohenberg derived an order parameter equation for the temperature and 
fluid velocity dynamics of the convection. This work has attracted much attention in the 
last decades and the order parameter equation has been interpreted as a model system 
for pattern formation [33]. For instance, Lega et al. modeled the spatio-temporal pattern 



formation in large aspect ratio lasers 
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371 ] by the SHE and Brazovskii 



38 



discussed 



the SHE in the context of the condensation of liquid crystals fluids. However, in spite of 
its broad applicability, the SHE fails to model some experimental results 33J and several 
previous studies extended the equation by adding some terms 39|, |40j. Hence the application 



of the SHE is not limited to the original physical problem and thus represents a generic model 
for the spatio-temporal dynamics of spatially extended systems j^. 

To gain the SHE from Eq. ([1]), we choose the temporal operator to T = d/dt and define 
h[V] = aV - bV 3 , g[V] = 0, f[V] = f x V, h < 0. Subsequently we obtain 

dV(x,t) 



dt 



h\V(x,t)]+fi [ dyK(x-y)V(y,t) 



(4) 
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In addition the kernel K{x) is symmetric and exhibits longer-ranged spatial interactions, 
i.e. K n — > for n > A. The subsequent application of Eq. (T5]) yields 

= aV(x,t) -bV 3 (x,t) + K Q f x V{x,t) 
By choosing b = K\\fx\jAK A > and re-scaling time and space according to t — > 



{AK A / Kl\f x \)t and x -> ^/2K 4 /K 2 x respectively we find the SHE [33] 

t] eV(x, t) - V\ X) t)-(l + ^X V(x, t) 



dt v ' ' v ' ' \ dx 2 



(5) 



with e = AK^aJ K^\fx\ — AK^Kq/ 1 K\ + 1. Hence, systems obeying the SHE neglect spatial 
interactions of order n > 6. 

Moreover, Eq. ([5]) maybe derived from an energy functional J-[V] by the generalized 
Ginzburg-Landau equation 

*L = -*aa ( 6) 

Since the present work aims to introduce nonlocal interactions to physical models in a gen- 
eralized way, we follow the approach of Ginzburg and Landau and attempt to model the 
spatio-temporal dynamics of spatially extended systems by a generalized energy functional. 
Taking into account the spatial interactions of elements in a system, the interaction energy 
in a spatial field between two spatial locations x and y at time t reads K(x—y)V(x, t)V(y, t). 
Here the term K(x — y) quantifies the interaction strength between both elements |4l|]. Con- 
sidering local interactions h[V] and their corresponding energy contributions W[V] addition- 
ally, then the field energy is the space integral over all spatial locations and all interactions 

T = l dxW[V(x,t\) 
Jn 

-~ [ dx [ dyK(x-y)V(x,t)V(y,t) (7) 
1 Jn Jn 

= I dxW[V{x,t)]-\v{x,t)^{-irKf^- 

J n n=0 

dxW[V{x,t)\ - \v{x,t) 
<> 2 

x ( K + K 2 \' + K A \' + ■ • 

ox 2 ax 4 



with — 5W/5V = h[V] and Q denotes the real axis. In the third line, we applied the identity 
([2]) . Then the specific choices applied in the previous paragraph yields the well-known energy 



functional of the SHE 
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J- = I dx^V 4 (x,t) 



V{x,t). 

(8) 



Corresponding to this derivation of the SHE from general physical considerations, we con- 
clude that the SHE represents a specific description of a general spatial system. This view 
angle on the equation is supported by its various motivations in physics mentioned above. 
In addition to the previous discussion, Eq. (J7j) is the Lyapunov functional of Eq. (j3J) as 

Then the question arises whether there exists an energy functional in the case of finite 
propagation speeds. Corresponding to the previous treatment, the first ansatz for a possible 
energy functional would contain the terms K (z)V (x , t)V (x — z,t — t(z)) with the distance 
z = x — y between two spatial locations and the propagation delay t(z) ^ taken from 
Eq. (0Q). Then applying the identity (j2J) in the time domain, we obtain 



K{z)V{x, t)V(x -z,t- t{z)) = 

■n—n \ / 



n=0 

We observe that the energy of a pair of elements at distance z, whose interaction is delayed 
by the finite propagation speed c, depends on the temporal derivatives of the elements. Sub- 
sequently the energy functional of the system T\V, dV/dt, . . .] would depend on the temporal 
derivatives of the field and the generalized Ginzburg-Landau equation (jBJ) is not applicable. 
Preliminary results on the different derivation of Eq. (J3J) including propagation delays by an- 
other variational principle seems promising. However, the discussion of these results would 
exceed the major aim of the present work and we refer the reader to forthcoming work. 
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C. Kuramoto-Sivashinsky equation 



At last we focus on the Kuramoto-Sivashinsky equation 33|,|43( which allows for the study 



44j, the 



of various phenomena in fluids and solids as e.g. the phase turbulence in fluids 
thermal diffusive instabilities of flame fronts 45j, the directional solidification in alloys [46] 



or the interface instability during the application of industrial beam cutting techniques [47 ]. 
Here, the terms in Eq. ([!]) are chosen as T = d/dt, g(V) = giV, f(V) = f 2 V 2 with gi, f 2 < 0. 
Now the kernel K(x) is non-symmetric and short-ranged, i.e. K n — > 0, n > 1 while the 
kernel L(x) is chosen symmetric. Hence the corresponding kernel moments obey L 2n+ i = 
and L(x) describes long-range interaction with L n — > 0,n > 4. Then the application of 
Eq. © leads to 

dV ft^ = {L gi-h )V{x,t) + L 2gi -^V{x,t) 
+ L 49l ^V( X ,t) + 2K lf2 V d ^ 

while h(V) has been chosen as h(V) = —K f 2 V 2 — h V, h > 0. Finally the re-scaling 
t — > (L4/L 2 \gi\)t and x — > \jL±j L 2 x yields 

= - vV - ^-V(x t) - ^V(x t) - 
dt V dx* l ' J dx± [ ' } dx 

(9) 



with rj = \gi\L + h and the additional condition f 2 = L 2 ^/ L 2 j L^j2K\\g\\. Equation (J9j) 
is called generalized Kuramoto-Sivashinsky equation |33| and r] = represents the original 
Kuramoto-Sivashinsky equation. Hence systems obeying the Kuramoto-Sivashinsky equa- 
tion involve symmetric spatial interactions up to order 0(a 5 ) and non-symmetric spatial 
interactions of order 0(a). It is interesting to note that this model includes both non- 
symmetric and symmetric spatial interactions while the reaction-diffusion equation and the 
Swift-Hohenberg equation involve symmetric spatial interactions only. 

D. Validity study 

After the previous examinations, the question arises which specific kernel functions may 
fulfill the appropriate conditions in the three models. To this end, we apply the Fourier 



S 



transform to the nonlinear function S(V(x)) in (jSJ) and obtain the spatial mode expansion 



with the sequence f n (k) = (—ik) n K n . In the following, we abbreviate /„ = f n (k) for conve- 
nience. As we will see later, f n (k) represents the contribution of interactions at order n to 
the spatial modes of the system. 

To obtain convergence of the series expansion in (j2J), the ratio of subsequent sequence terms 
\f n \ has to be smaller than unity. This means this condition reads |/ n +2|/|/n| < 1 for n even 
and |/n+i|/|/n| < 1 f° r & h n f° r symmetric and non-symmetric kernels, respectively. 
To be more specific let us focus on a specific symmetric kernel function. Since the kernel rep- 
resents the probability density function of the spatial interactions between two locations, the 
central limit theorem supports the choice of the Gauss function K{x) = exp(— x 2 /2a 2 ) / \f2jio 
with the spatial range constant a. Then the previous convergence condition leads to 
a 2 < (n + 2)/k 2 ,n = 2,4,6,... . If there is a maximum spatial frequency \k m \, i.e. the 
Fourier transform S'( | A; | ) is negligible for \k\ > \k m \, then the spatial interaction range a is 
limited from above by y/n + 2 / 1 k m \ . Subsequently, the spatial interaction range a is delim- 
ited to the interval < a 2 < l/k^ which guarantees the Taylor expansion convergence in 
d2J). We point out that the limitation to spatial modes is present in various physical systems. 
For instance, spatial systems near a phase transition exhibit a hierarchy of time scales of 
spatial modes with corresponding prominent spatial frequencies. In this case the time evo- 
lution of the Fourier transforms S(k,t) at spatial frequencies far from the prominent ones 
are negligible and the prominent frequencies define the maximum spatial frequency \k m \. 
However, if there is no maximum spatial frequency, i.e. \k m \ — > oo, it is a — > and the con- 
vergence condition yields K(x) —>■ 5(x). Hence for an unlimited range of spatial frequencies, 
the expansion (|2j) is valid for local interactions only and K n — > for n > 0. 
Now let us interpret the latter results. It turns out that the convergence condition for 
Gaussian kernels leads to the dependence of the maximum spatial range on the interaction 
order n. For instance if the expansion in Eq. (J2J) approximates local interactions but no 
higher orders of nonlocality (cf. the previous discussion of the reaction-diffusion equation), 
it is \K n+2 /K n \ <C 1 for n — 2,4, . . . and, subsequently, a 2 < (2 + 2)/k 2 = 4/k 2 . Hence, 
a maximum spatial frequency \k m \ yields the validity condition < o 2 < In case of 




oo 



n=0 



(10) 
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spatial interactions at higher orders of nonlocality (cf. the previous discussion of the Swift- 

Hohenberg and the Kuramoto-Sivashinsky equation), it is \K n+ 2/K n \ <C 1 for n — 4, 6, 

This means that the validity condition now reads < er 2 < (4 + 1)jk 2 m = 6/k^. Thus the 
spatial interaction range a may be larger than in the diffusion model. 
Since the Kuramoto-Sivashinsky equation involves non-symmetric interactions, let us also 
examine briefly the kernel function K{x) = 5{x — xq) with the non-symmetric shift xq ^ 0. 
We find K n = x^/nl and the convergence condition reads x < (n + l)/\k m \. Subsequently, 
the Kuramoto-Sivashinsky equation is an approximation of ([I]) if Xq < 2/\k m \ and if a max- 
imum spatial frequency exists. 



IV. LINEAR ANALYSIS 



Now let us examine the stationary state and its linear behavior for a finite propagation 
speed c. In case of the stationary state V(x,t) = V , Eq. ([1]) recasts to 

T(0)V = h(V ) + K K f(V ) + K L g(V ), (11) 

/oo poo 
K(x)dx, k l = I L(x)dx. 
-oo J — oo 

with k k = J K(x)dx, k l = J L(x)dx. Considering linear deviations u(x,t) = V(x, t) — V ~ 
e xt+ikx : Eq. (HJ reads 

/oo 
dzM(z)e- xlzl/c e- ikz (12) 
-oo 

with M(z) = SfK(z) + s g L(z) and Sh = Sh/5V, Sf = Sf /SV, s g = Sg/5V computed at V = 
Vq, while 5/SV denotes the functional derivative. In the following, the propagation delay c is 
assumed being large but finite and the approximation exp(— A|z|/c) ~ 1 — A|^|/c + A 2 (|^|/c) 2 
holds. For the characteristic spatial scale of the system a, we introduce the characteristic 
propagation delay r = a/c which represents the delay time caused by the finite propagation 
speed. Thus the previous approximation holds for small propagation delays r <^ 1/A. 
Subsequently expanding the exponential exp(— ikz) into a power series we obtain 

-. oo \ 2 00 

T(X) + - c J2(-lk) n P n - ^ ^("^) n Qn 
n=0 n=0 

oo 

= S h + J2(-* k T M n (13) 

ra=0 
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with the kernel moments M n defined previously and P n = J n dzM(z)\z\z n /n\, Q n = 
J n dzM(z)\z\ 2 z n jn\. Preliminary computations of M n , P n and Q n for Gaussian kernels yield 
M n ~ 0(a n ), P n ~ 0(a n+1 ) and Q n ~ 0(<7 n + 2 ) 

Equation ([TBI defines the stability condition for general spatial interactions and short prop- 
agation delays. Additional more detailed studies show that large propagation speeds c with 
1/c 2 ~ and vanishing odd kernel moments M 2n +i = 0, Pm+\ =, Qin+\ = lead to real 
Lyapunov exponents and thus may yield stationary bifurcations. In contrast oscillatory bi- 
furcations may emerge for non-vanishing odd moments only. By recalling the definitions of 
these moments, we find that symmetric and non-symmetric spatial interactions may yield 
stationary and non-stationary bifurcations, respectively. In case of smaller propagation 
speeds 1/c 2 7^ 0, this classification does not hold anymore and oscillatory bifurcations may 
occur for both symmetric and non-symmetric kernels. 

A. The reaction-diffusion equation 

To illustrate the latter results, first let us study a specific react ion- diffusion equation, 
namely the homogeneous diffusion equation also known as Ficks Second Law. According 
to the previous discussion, Eq. (j3J) represents the diffusion equation with h = yielding 
h(V) = —fiV. Further Eq. ( TTTT) yields the stationary solution Vo = and the parameters 
Sf = fi, Sh = —fx, while the spatial interactions are given by K(z) 7^ and L(z) = 0, i.e. 
M{z)=f 1 K{z). 

According to the previous discussion of fflUj) . Eq. (j2J) represents a spatial mode expansion 
in orders of the characteristic spatial scale a. That is, in the case of the diffusion equation, 
the order expansion is truncated at 0(a 3 ) yielding M — fx ~ 0(1), M 2 = /icx 2 /2 ~ 0(a 2 ) 
while all odd and higher kernel moments Mm vanish with n > 2. Further we find Po = 
/iv^/v 7 ^, Pi = 0, P 2 = /iV2<7 3 /0r with P n w V n > 3 and Q = fca 2 , Q n ^0\/n>0. 
Here it is f\ = 2D a 2 and D represents the diffusion constant. Subsequently in the case of 
c — >• 00 Eq. (flBl yields X(k) ~ —Dk 2 which is the well known dispersion relation of the 
diffusion equation. Figure [1] compares A and A#£> and the corresponding resulting field 
activity for two values of a. We observe that the IDE model is a good approximation to 
the diffusion equation for a = 0.1, while the the approximation is much worse for a = 2. 
This finding confirms our previous results on the convergence condition. Here the maximum 
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field u(x,t) 



a=2.0 



a=0.1 



wave number k 




-5 5 

space x 



FIG. 1: (Color online) The field of the diffusion equation and the Lyapunov exponent for different 
spatial ranges a. The dashed lines represent the results from the diffusion PDE, while the solid lines 
denotes the results from the IDE model. In case of a = 0.1, the reconstructed field u(x, t) obtained 
from the IDE can not be distinguished from the field of the original diffusion equation. In addition 
the colors in the left panel denote subsequent times t = (black, very thick line), t = 3 (red, thick 
line) and t = 6 (green, thin line). That is the field u(x,t) flattens at increasing time. 

value of \k\ is 2/ a according to the convergence condition a < 2/k derived in the previous 
paragraphs. 

Now let us turn to finite propagation speeds c < oo. Considering the corresponding 
orders of M n , P n and Q n , Eq. f[T3j) reads 



A + -(P - P2k 2 ) - %Qv = s h + M - M 2 k 2 . 



(14) 



For small propagation delays (A7 



>^RD{k) 



we find the Lyapunov exponent 
—Dk 2 



1 + ar(l - a 2 k 2 ) 



(15) 



with a = fiy/2/ir. By virtue of the wave number limit \k\ < 2/ a, there is a critical 
propagation delay r th = l/3a and a corresponding critical propagation speed c t h = Saa. 
Then small delays r < r t h retain the field stability, while large delays r > r t h yield linear 
unstable modes k with A^d(^) > 0. Figure [2] shows the corresponding stability diagram of 
the system. 

Moreover, Eq. (|15|) may be interpreted as if it originates from the modified diffusion 
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wave number k 



FIG. 2: Illustrated stability diagram of the extended diffusion equation. The vertical axis k = 1 
represents the asymptote of the instability regime. 



equation 



du(x,t) aa f 2 d 2 ^ du(x,t) 
c V dx 2 



D 



d 2 u(x, t) 



(16) 



dt c \ dx 2 J dt dx 2 

This equation is an extension of the standrad diffusion equations considering finite propa- 
gation speeds. It applies in media whose propagation delay on the typical length scale of 
the system is not negligible and, for instance, which show delayed spread of space-time ac- 
tivity. We mention results from fast hot pulses in plasma [11], whose investigation indicates 
nonlocal effects in the reaction-diffusion mechanism. Moreover, there is the debate on the 
non- Fourier heat conduction in nonhomogeneous materials which show such delayed tempo- 

nnn 

ral activity at measurement points [8|, [9J, |48[ . In order to compare our model to the results 
found in experiments Fig. [3] shows the temporal activity at a spatial point. We observe 
a prominent delay of the activity at large delays, while r — > yields the Fourier case. This 
finding is irrespective of the spatial location used and coincides qualitatively to experimen 
tal findings. Moreover, the derived model (fl6|) contrasts to the Cattaneo equation 
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491, 



which represents a telegraphic equation and extends the parabolic diffusion equation to a 
hyperbolic partial differential equation by an additional second temporal derivative and thus 




time t 



FIG. 3: Temporal activity of the modified diffusion equation at spatial location x = 8 for two 
different propagation delays, i.e. propagation speeds. The solutions are computed analytically 



with initial condition u(k,t = 0) = exp(-/c 2 /2 • 0.8 2 )/v / 2^ra8 f 



13 



involves finite propagation speed. Originally the Cattaneo extension has been proposed ac- 
cording to heuristic arguments while our model is derived from the explicit treatment of 
propagation delay. 

At last, let us briefly focus on the case of larger propagation delays (Ar) 3 — > yielding the 
implicit equation —\ 2 t 2c 2D j o 2 + A(l + aDr(l — o 2 k 2 )) — Dk 2 = 0. Simple calculus reveals 
Re(X) > for some k and larger delays destabilize the field similar to the previous case of 
small delays. 



B. The Swift-Hohenberg equation 

Now let us turn to the Swift-Hohenberg model which involves long-range spatial inter- 
actions. It has been studied as a generic model for spatial pattern formation in various 



physical systems [36|, |38|, |40|. By virtue of its broad applicability, this paragraph aims to 



extend it by finite propagation delay effects. The stationary state is given by Eq. ( TTTT) and 



we find one or three solutions. The single solution Vo = leads to the parameters sj = —\fi\ 
and Sh = a. Further the spatial interactions are given by K(z) ^ and L(z) = yielding 
M(z)=-\f 1 \K(z). 

According to the discussions in section UlIDl Eq. (J2]) breaks off at order 0(a 5 ) and we obtain 
the expansion terms 

M = -IM, M 2 = -\h\a 2 /2, M 4 = -|/i|a 4 /8 

2 . ,,, /2 3 D _ V2 * 



P = -\fl\\ -°, P 2 = -i/lU/-^, P 4 = H/li^tf 
7T V VT 3V7T 

Qo = -l/ik 2 ,g 2 = -|A|^ 4 

while all odd and higher terms M n , P n , Q n vanish. Subsequently Eq. ( fT3l) yields for c — > oo 

Aoo = ^(e-l + 2/ 2 -/ 4 ) 

with the scaled wave number I = akj \/2 and e taken from section [Til Bl Here Aqo represents 
the Lyapunov exponent of the linearized original Swift-Hohenberg equation (jHJ) and thus 
defines the linear stability of the underlying spatial system. To be more detailed, for e < the 
system is linear stable with < for all I, while e > yields Aoo > for \—yfe < I < 1+y/e. 
In other words, a Turing instability emerges for e > at scaled wave numbers |/| ~ 1. 
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Now let us turn to the case c < oo. For small propagation delays (At) 2 — > we find 

XsH = l-a(l-2F + 4/ 4 /3)r ^ 

with a = |/i|a/2/7t. Here As# represents the Lyapunov exponent of the Swift-Hohenberg 
equation subject to large propagation speeds. It is important to mention that the validity 
criteria < a 2 < derived in section Ull Dl delimits the scaled wave number to —y/3 < 

I < y/3. Taking into account this constraint, closer examinations of Eq. (|17p reveal the 
critical propagation delay r t h = I /7a and the corresponding critical propagation speed 
c t h = 7aa. For r < r^, the denominator of \$h is positive for all valid wave numbers and 
thus Aoo defines the stability of the system. Hence the propagation delay does not affect the 
stability of the system. 

In contrast large propagation delays r > yield a negative denominator for a half-band 
of wave numbers |/| > l t h and thus change the sign of the Lyapunov exponent. Here the 
critical wave number l t h represents the root of the denominator in Eq. (IITI) . Figure H] shows 
the stability diagram of Eq. (IT7|) . We observe the unstable half-band for r > for all e, 
while t > 1/a > r t h yields an additional unstable band of wave numbers around I = 0. For 
e < the stable band of spatial modes vanish for large propagation delays r > 4/a and the 
system is totally destabilized, while for e > increased propagation delays stabilize the band 
of unstable Turing-modes. This stabilization by propagation delays represents a novel effect 
near Turing instabilities. In addition recall that the previous treatment around l t h reflects 
|A| — > oo and thus conflicts with the previous assumption (Ar) 2 rs 0. Subsequentially higher 
polynomial orders in A are necessary for closer examinations. 

Finally Eq. (117j) may be interpreted as if it originates from the extended Swift-Hohenberg 
equation 

dV(x,t) aa A | 2 d 2 4^^ dV(x,t) 



dt c \ dx 2 3 dx 4 J dt 

Q2 



= eV(x, t) - V\x, t)-\l + w \ V(x, t) (18) 

with the re-scaled time t — > 2i/|/i|. This equation represents an extension of the well-knwon 
Swift-Hohenberg equation involving a finite propagation speed. In the case of a infinite 
propagation speed c — > oo the additional second term on the left hand side of Eq. ( fT8l) 
vanishes and the original model equation is re-gained. Further the finite propagation speed 
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£ < 




scaled wavenumber 



FIG. 4: Illustrated stability diagram of the extended Swift-Hohenberg equation involving propa- 
gation delays. For e < (top panel) the underlying system is stable for r < l/a(l — 2l 2 + 4/ 4 /3) 
and unstable otherwise. For e > (bottom panel) the stability is the same as for e < besides 
the wave number band e_ < 1 < e +) s± = ±yfe. For short delays r < 1/q(1 - 2l 2 + 4Z 4 /3) the 
underlying system is linear unstable in this wave number band, while the modes in this band are 
stable otherwise. In both latter cases r < r-th guarantees the linear stability of the underlying 
system. 

leads to the coupling term between spatial and temporal derivatives in Eq. ffl~8l) reflecting 
the space-dependent delay in Eq. (CQ). 



C. The Kuramoto-Sivashinsky equation 

Finally let us investigate the Kuramoto-Sivashinsky equation by following the same steps 
as in the previous examples. Now Eq. (jl~3l) reads for small propagation delays (At) 2 w 

A + -(P - k 2 P 2 ) = s h + M - ikM 1 - k 2 M 2 + k A Q 2 . 
c 

with M = s f + s gi M x = SfX Q , M 2 = s g a 2 /2, M 4 = s 9 a 4 /8, P = s f x + s g ^2a/^, P 2 = 
s g \^2a 3 / y/ft and M n , P m — >• for n > 4, m > 2. The stationary solution is taken from ffTTj) 
and reads Vq = yielding s/ = 0, s g — g\ and the Lyapunov exponent 

-n _)_ 72 _ 74 

A ^ = 1 Vl 1^ ( 19 ) 

1 — ra(l — 4r) 
16 



with a = \g\\ y2/V> l = &k/2 and Axs = A/2 1 <7i | . We point out that the re-scalings k — > / 
and A — > Xks are identical to the re-scalings of space and time in section IIII CI 
Now let us examine ([T§j) in some more detail. For a vanishing propagation delay r = Xks 
becomes the Lyapunov exponent well-known for the original Kuramoto-Sivashinsky model. 
This means the stationary solution is linear unstable at wave numbers — 1 < I < 1 while it 
is linear stable otherwise. In case of r ^ 0, the propagation delay changes the Lyapunov 
exponent and may destabilize the system. It turns out that there is a critical propagation 
delay r t h = 1/at and a corresponding a critical propagation speed c t h = cto. For small 
propagation delays r < r t h, the Lyapunov exponents are modified but the system stability 
is retained, see Fig. In contrast large propagation delays r > r t h lead to a sign inversion 
of the Lyapunov exponent for — l c < I < l c with I 2 = 1 — 1/ err. Figure [6] summarizes the 
stability conditions in a diagram. However r — * r t h yields A — > oo for / — > l c and thus does 
not fulfill the previous condition (At) 2 w 0. 

Finally, Eq. (fl9l) may be interpreted as if it originates from the PDE 

dV(x,t) aa f l ^d 2 \dV{x,t) 



dt c \ dx 2 J dt 

d 2 d 4 dV 

_^__ ni , t) __ nx , f) _y_ (20 ) 

This equation represents an extension of the Kuramoto-Sivashinsky equation now involv- 
ing large propagation speeds c. We observe that c —>■ oo yields the original Kuramoto- 
Sivashinsky equation, while the finite propagation speed leads to the additional second term 
on the left hand side. 




FIG. 5: The Lyapunov exponent of the Kuramoto-Sivashinsky equation subject to the scaled wave 
number /. The different line styles encode the different propagation delays r = (dotted line), 
r = 2.0 (dashed line) and r = 2.499 (solid line). It is T th = 2.5. 
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unstable 

-0.5 0.5 1 

scaled wave number I 



FIG. 6: The illustrated stability diagram of the extended Kuramoto-Sivashinsky equation subject. 
The center stability region does not exceed the interval [—0.5, 0.5] on the horizontal axis. 

V. CONCLUSION 



The present work has introduced a nonlocal integral-differential equation model which 
generalizes three important partial differential equations. Moreover we have shown how 
to consider large propagation speeds in these models and have quantified the propagation 
delay as the fraction of spatial interaction range and propagation speed. The subsequent 
linear stability analysis turns out to be dependent on this quantity and reveals critical 
propagation delays and corresponding propagation speeds. We find that large propagation 
speeds slow down the activity spread in diffusion systems and thus allows for the expla- 
nation of non-Fourier behavior in non-homogeneous systems. Further finite propagation 
speeds may destabilize the dynamics of Swift-Hohenberg equation in the stable regime of 
the original equation, while it may stabilizes an occurring Turing instability. Finally finite 
propagation speeds may stabilize the dynamics of the Kuramoto-Sivashinsky equation. For 
all three models novel partial differential equations have been formulated which incorporate 
the propagation delay effects. We are convinced that propagation delays play an important 
role in media showing ultra-fast phenomena and that the obtained results shall allow for 
their modeling. 
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